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Abstract 

We perform a self-consistent relativistic RPA calculation for the 
isobaric analogue and Gamow- Teller resonances based on relativistic 
mean field theory results for the ground states of ^^Ca, ^"Zr and ^'^'^Pb. 
We use the parameter set NLl for the a, uj and p mesons, and exper- 
imental values for the pion and nucleon. An extra parameter, related 
to the intensity of the contact term in the pion-exchange interaction, 
is crucial to reproduce the latter resonances. 
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In recent years the relativistic mean field theory (RMFT) has been suc- 
cessfully applied to account for both: i) the bulk features of nuclear matter 
(saturation, equation of state, etc.) and ii) the ground state properties of 
finite nuclei, including unstable ones up to the nucleon drip lines [0]. (For 
reviews, see refs. 0, |^, |^.) A nice feature of the model is that the 
same restricted set of parameters can be used for all these systems. Even 
though some excited nuclear states, including certain giant resonances, have 
also been calculated in the relativistic random phase approximation (RRPA) 
based on RMFT, this has never been done for charge-exchange excitations. 
The aim of the present work is to explore to what extent this relativistic 
model, whose parameters have been fitted to ground state properties, can 
describe the charge-exchange collective states, such as the isobaric analogue 
state (IAS) and the Gamow- Teller resonance (GTR). 

Both IAS and GTR have been extensively calculated in the context of 
traditional, i.e., nonrelativistic, nuclear structure theory, usually in the ran- 
dom phase approximation (RPA) or extensions of it, making use of, not 
only completely phenomenological interactions, but also more realistic ones 
based on meson exchange. While the IAS is easily reproduced, the GTR is 
more problematic, especially with respect to its strength that turns out to 
be quenched to ~ 60% of the Ikeda sum-rule value. (For a review of GTR, 
see ref. [0.) Two mechanisms have been proposed for the absorption of the 
unseen strength: i) the excitation of the A resonance, which has the same 
quantum numbers as the GTR, or ii) the influence of 2p-2h and higher exci- 
tations, not included in RPA. We shall not dwell on such finer points in this 
paper, since this would require, anyway, a more sophisticated treatment than 
the simple version of RRPA we are using, but rather will concentrate on the 
quality of the reproduction of excitation energy and strength as compared 
with similar nonrelativistic calculations [p!2| , p^ . 

Relativistic mean field theory (RMFT) always takes as a starting point 
some version of an effective quantum field theory describing a system of inter- 
acting hadrons, coupled also to the electromagnetic (or to the electroweak) 
field, generally referred to as quantum hadrodynamics (QHD) ^, ^, |]. 
Historically, QHD was introduced as a renormalizable quantum field the- 
ory, which severely limited the number of admissible terms in its Lagrangian 
[0]. However, despite its many successes, it had difficulties with its renor- 
malization program. Thus, in the last years, it began to be formulated as 
an effective, nonrenormalizable quantum field theory representing the low- 
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energy limit of the fundamental theory of strong interactions, namely, quan- 
tum chromodynamics (QCD). Consequently, all terms compatible with the 
symmetries of QCD were now admissible, and one had to find appropriate 
expansion parameters and a criterion of naturalness to select the most im- 
portant ones. Even thus, the number of terms increased significantly. Their 
coupling constants were adjusted to a judicious selection of nuclear properties 
that should be well reproduced at the mean field level and, while this was 
successful in the sense that they turned out to be natural, on the other hand 
there are clear indications that their full set is underdetermined |^. Since 
our main purpose here is not to test the foundations of QHD, but rather to 
assess its ability to reproduce some new nuclear properties, and we want to 
keep the number of free parameters to a minimum, we will include only the 
most important terms, and work with the following nonrenormalizable QHD 
Lagrangian density: 

1 111- 

- \f,,F^'' - e^7^1±^^A^ (1) 

in standard notation {% = c = 1), where ip is the spin-|, isospin-| nucleon 
field, cr, o;^, tt and denote the meson fields, and A^^ is the electromagnetic 
field 0. The field tensors for the vector particles read 

To reproduce the experimental values of the nuclear incompressibility and 
surface diffuseness in RMFT, it has been found important [§] to allow the 
sigma meson to self- interact. This is done through the cubic and quartic 
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terms, in cr, appearing in eq. (|1]). All couplings are taken to be of the direct, 
nonderivative kind, except in the pion-nucleon case, where both the phe- 
nomenology and theoretical reasons (related to chiral symmetry) definitively 
favor the pseudovector coupling. It is true that phenomenology would also 
recommend that a certain amount of derivative coupling should be mixed in 
for the vector mesons to the nucleon Yet, for simplicity, we followed the 
general trend in RMFT calculations [0, 0, ^ and did not include such terms 
in eq. (|I]). For the omega meson, this can be justified since that term would 



be rather small |I0[, but for the rho meson, while unimportant for the nu- 
clear ground state, the corresponding term would give a sizable contribution 
to excitations of the Gamow- Teller type. Finally, the photon is minimally 
coupled to the proton to take care of the important Coulomb interaction 
inside nuclei, but not to the charged mesons since this would be of much 
less relevance for the properties of interest here and is hence ignored in this 
model Lagrangian. 

The nucleon mass Q, M, the proton charge, e, the pion mass, m^, and 
coupling constant, are fixed at their experimental values. Thus, the 
masses of the remaining mesons included in eq. mg-, and rrip, their 
coupling constants, g^, Qu) and Qp, and the self-interaction strengths for the 
sigma meson, g2 and g^, give a total of only 8 free parameters for this model. 

RMFT results from an approximation scheme to solve the Euler-Lagrange 
equations derived from £, which consists in taking advantage of the high 
values of the densities inside nuclei to replace the meson and photon fields 
by classical fields satisfying these equations with the source terms replaced 
by their expectation values in the nuclear ground state of interest. This is 
equivalent to compute the meson self-energies in the Hartree approximation 
disregarding the negative-energy, i.e., antiparticle, states, which is sometimes 
referred to as the no-sea approximation. 

The RMFT equations simplify considerably if one takes advantage of 
known symmetries of the nuclear ground state. Firstly, as it always has 
a definite parity, the pion, due to its pseudoscalar nature, completely dis- 
appears from these equations. Secondly, the definite charge of the nucleus 
eliminates all but the third component of the isovector mesons. Also, we 
are interested here only in spherical nuclei, and their rotational invariance 
implies that the spatial components of the currents in the source terms van- 



^The proton mass is used here. 



4 



ish in RMFT. So, the corresponding components of the vector fields can be 
ignored. Finally, we are looking only for stationary solutions in the nuclear 
rest frame. Thus, the meson and photon fields are taken as time-independent 
and the nucleon field has merely time-dependent phase factors, i.e., 

i'i^f^) = H (^cMc{r) expi-iEat), (3) 

a 

where the no-sea approximation has been made, and Oq, is the annihilation 
operator in the particle state with (positive) energy E^- With these simpli- 
fications, the RMFT equations become 

Ha- V + /5(M + K(r)) + K(r)]Wa = EMa, (4) 

and 



(-V^ + m^2)cr = -gaps{r) - g2(T'^ - g-i(T^, 
(-W^ + mj)uj'' = g^p,{r), (5) 
(_V2 + m/)p° = gppsir), 
-VM° = epp(r), 

where 

Vs = g.a, V,=g^uj^ + g^Tspl + e^^^A'>, (6) 
are the scalar and vector potentials, and 

Ps = (^Pi^) , Pv = (V^tV) , P3 = (^7°r3^) , Pp = (^7°^^V^) , (7) 

the scalar, vector, isovector (third component) and proton densities. The 
expectation values are taken in the nuclear ground state, constructed by 
putting Z protons and neutrons, {Z + N = A), in the lowest particle 
states with equal to +1 and —1, respectively, obtained by solving the 
Dirac equation (^). This yields the following results for the densities 

A _ A _ 

Ps = ^UJAa, Pv=^l^al^Ua, 

P3 = E^^T^sWa, Pp= E^«7°^^W„. (8) 

a=l Q=l ^ 
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The set of eqs. (H)-(|D should be solved self-consistently. 

Alternatively, the RMFT equations could have been obtained by extrem- 
ization of an energy functional built from eq. and involving only valence 
nucleons and classical photon and meson fields. This density functional ap- 
proach is sometimes preferred when working strictly within the context of 
an effective field theory, since it becomes clear in this case that much of the 
effects of vacuum dynamics and of certain many-body correlations beyond 
the Hartree level are absorbed in the renormalized coupling constants in the 
Lagrangian, when fitted to experimental data [||, 

The RRPA equations have been derived several times in the literature and 



applied, both to nuclear matter |]T6|, and to finite nuclei JT^^lMIl- -^^^ '^^^ 
purpose here, they can be obtained from the equations of motion formalism 
of Rowe [^. One introduces charge-exchange excitation operators of the 



form ra, El 



pn np 



where p and p {n and n) label unoccupied and occupied proton (neutron) 
positive-energy, single-particle RMFT states. J is the total angular momen- 
tum, and A runs from 1 to vr + z/, with vr and u being the number of pn and 
np excitations, i.e., those having J^pn \^pn\ ~ ^np I'^npl equal to +1 and 
— 1, respectively. Following the standard procedure one arrives at the RRPA 
equations 

(^^t ) (y^A ) =^Ja(^,a ). (10) 
The submatrices are given in terms of the residual interaction V by 

^in,p'n' = (Ep - En)6pp:6nn' + {{pn-^)j\V\{p'n'-^)j), 

= {E^- Ep)5nn'5pp+ {{np-')j\V\in'p'-')j), (11) 
Binyp' = {-y^'''{ipn-')jMAn'r')j,-M,\V\). 

The excitation energies Ej^^ (A+ = 1, ■ ■ - vr) in the odd-odd {N — 1, Z + 1) 
nucleus (measured from the target ground state) are the highest vr eigenvalues 
u j\ in (|13). In the same way the excitation energies -Exx_ (A- = 1, ■ ■ ■ z/) in 
the (A^ + 1, Z — 1) nucleus are the remaining eigenvalues Ujx taken with 
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opposite sign. The corresponding total transition strengths are: 



S1 



Si 



E 

A+ 

E 



J2Xpn^*{p\\wj\\n) + J2Ynp''^*{p\\u'j\\n) 



JA_|.* 



pn 



np 



^Ynp {p\\wj\\n) + ^Xp^ {p\\wj\\n) 



JA_ 



np 



pn 



where Wq = 1 and Wi 



= (J. They fulfill the well known sum rules: 



(12) 



(13) 



For a self-consistent calculation, V must be obtained from the same La- 
grangian, given in eq. (P, which was used for the mean field part. Also, 
to be compatible with RMFT, one must consider only the direct terms in 
eq. (11^)? leading, in fact, to a ring approximation. Hence, only the isovector 
mesons contribute to isovector excitations like IAS and GTR. Furthermore, 
when the instantaneous approximation is made, the interaction, V = Vj^ + Vp, 
reads 



v;(i,2) 



— Ti ■ T2 (cTi ■ Vi cr2 ■ V2)y{m^, r). 



V^p(l,2) = Qp^Ti- T2{1 - Oil- OL2)Y{mp,r), 



(14) 



with r = |ri — r2| and Y{m,r) = exp(— mr)/(47rr). 

We performed our numerical calculations taking '^'^Ca, ^^Zr and ^''^Pb as 
the target nuclei. To get a discretized set of single-particle RMFT states, we 
solved eqs. (^-(H) by expanding the different fields in truncated harmonic os- 
cillator bases, following the method of P. Ring and collaborators as described 
in ref. 0, where the details can be found. For the RMFT parameters we 
chose the set NLl, which is widely used for nuclear ground state calculations 
[0, 1^, ^. The RRPA equations were solved for J'" = 0"*" and 1"*" states in a 
model-space including only OhQ and 2hQ excitations. Furthermore, since the 
continuum is not well represented by an expansion in a harmonic oscillator 
basis, we tried to avoid it as much as possible. To this end, we included 
only those RMFT single-particle states that are bound at least for neutrons. 
(Of course, for consistency, we had to include the corresponding states for 
protons, even when unbound.) That this model-space is not too restricted is 
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attested by the fact that, in all the cases considered here, the sum rule (0) 
was obeyed at the level of 99.5% (92.8% ) or better for the IAS (GTR). 

Our results for the IAS of ^osp]-, gjyg excitation energy (always mea- 
sured from the target ground state) of 18.6 MeV and a strength of 99% of 
the sum-rule prediction of (A^ — Z), which compare well with the experimen- 
tal values of 18.8 MeV and ~ 100%, as quoted in [^. These results are of 
the same quality as those of Colo et al. 0], who perform a nonrelativistic 
calculation that is also self-consistent, in the sense that they use the same 
interaction (Skyrme force) for the excitations as for the target ground state 
and single-particle basis. 

Since the GTR's are rather broad, having experimental widths of ~ 
4 MeV, the theoretical results we report here correspond to a smoothed 
strength function, constructed by replacing the spikes in the RRPA strength 
distribution by Lorentzian peaks of conveniently chosen width. The excita- 
tion energy and strength of the GTR are then obtained by fitting a Lorentzian 
to the main peak in this strength function. When the interaction (0) is used 
we get that the excitation energy and the strength of the GTR in ^osp]^ ^^.^^ 
respectively, ~ 6 MeV and 56% of the Ikeda sum-rule, while the correspond- 
ing experimental values are 19.2 MeV and 60 — 70% [0. Therefore, as it 
stands, the relativistic calculation would underestimate the GTR by more 
than 10 MeV, putting it much lower than the IAS, in blatant opposition to 
experiment. The nonrelativistic calculation of Colo et al. does much better, 
though it overestimates the excitation energy by 2 — 3 MeV. 

It is still possible, however, to improve the GTR results in our calculation, 
without losing the self-consistency between the excitation and the mean field 
parts. In fact, when dealing with the pion contribution, we might follow 
the general practice and eliminate the contact term, which comes from the 
derivative (pseudovector) coupling and is thought to be suppressed by the 
short-range correlations between nucleons. (As such it should disappear when 
a form factor is used to account for finite nucleon and meson sizes.) More, 
to reproduce the energetics of the GTR within a nonrelativistic calculation 
based on tt- and p- exchanges, Krewald et al. |T^ were forced to add a 
zero-range Landau-Migdal interaction of adjustable strength. In our case 
this would correspond to gauge the contact term in the bare 7r-exchange 
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interaction, given by (p!4D, or equivalently to add the term 



2) = g' Ti ■ T2 o-i ■ 0-2 (5(ri - r^), (15) 



which for g' = 1/3 completely cancels the contact part in ([T^)- This term was 
found to play an important role in the RRPA calculation of unnatural parity 
states in ^^O made by Blunden and McCorquodale |T^, who used the purely 
phenomenological value g' = 0.7, in analogy to nonrelativistic calculations. 
Since we are neglecting exchange contributions, has no effect on the IAS, 
but is of overwhelming importance for the GTR, as illustrated in Fig. |I] for 
the case of ^°^Pb. One sees that the theoretically justified coupling g' = 1/3 
leads to better results than the bare interaction {g' = 0). Yet, the agreement 
with experiment is only achieved with g' = 0.7. Thus, we found it necessary 
to incorporate 5Kr, with g' = 0.7, into the residual interaction. 

The results obtained in our calculation, with this choice of residual in- 
teraction, for the IAS and GTR, excited from the closed-shell nuclei ^^Ca, 
^°Zr and ^''^Pb, are summarized in Table |1]. One can see that they are of 



the same quality as those obtained in similar nonrelativistic calculations [14 



The agreement with experiment is very good for the excitation energies of 
both resonances, but less so for the strengths, specially for the GTR. This 
was already expected, since our model does not include either 2p-2h or A-h 
excitations, and therefore is unable to reproduce the quenching of the GTR. 
For similar reasons, and also because we do not treat the continuum prop- 
erly, we get very little strength at higher excitation energies beyond the GTR 
peak. In the case of ^°Zr, for instance, our calculated strength function is 
decaying well below 0.45 MeV~^ already for ~ 25 MeV excitation, while the 
experimental one p9| stays at about this value up to ~ 60 MeV. 

We conclude that both the IAS and the GTR excited from closed shell 
nuclei can be well reproduced in the context of quantum hadrodynamics, 
with the ground state parametrization for the a, uj and p mesons, and ne- 
glecting the exchange contributions. For the GTR, however, the pion with 
pseudovector coupling to the nucleon also must be included, and it is essential 
to use for the strength g' of the contact term a value close to that used 
in similar nonrelativistic calculations. An alternative approach, to account 
for the energies of the GTR's, might be the mixed coupling of the rho meson 
to the nucleon. Yet, such a parametrization for ground state properties is 
available only in the Hartree-Fock approximation . 
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Table 1: RRPA results for excitation energies and strengths of charge- 
exchange resonances calculated with isovector interaction including a 
Landau-Migdal-type contact term of intensity g' = 0.7, as described in the 
text. Experimental values are given for comparison. 



Parent 




RRPA 


Experiment 


and 


Energy 


Strength 


Energy Strength 


resonance 


[MeV] 


[% of sum rule] 


[MeV] [% of sum rule] 



208pbg 

IAS 18.6 99 18.83 ±0.02 -100 

GTR 18.9 80 19.2 ±0.2 60 - 70 

90Zif| 

IAS 11.9 101.3 12.0 ±0.2 -100 

GTR 16.0 80 15.6 ±0.3 28 

4«Ca[| 

IAS 7.48 100.4 7.17 ~ 100 

GTR 11.0 82 ~ 10.5 35 



"Experimental values taken from Refs. 



''Experimental values taken from Refs. |28, 2£]. 
'^Experimental values taken from Ref. 



271, 
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Figure 1: Theoretical Gamow- Teller strength distribution for the parent 
nucleus ^osp]^ obtained with an isovector interaction including a Landau- 
Migdal-type contact term of strength g' = 0, 1/3 and 0.7. The spikes give the 
raw RRPA results and the continuous curve, the strength function smoothed 
out as explained in the text. The strength function for the resonance peak 
extracted from experiment |^ is drawn in dotted line. 
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